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OPTIMUM PREDICTION FOR A 


TORQUE SATURATING SERVO 


SUMMARY 


This study was made to examine the use of a particular type of nonlinear 
control for a torque saturating second order servomechanism. The control 
involves mechanization with diode limiters and operational amplifiers of a 
torque switching curve in the phase plane of error and error rate, with a 
region of linear operation about the curve. The switching curve is shaped to 
give prediction for optimum time response and essentially no overshoot for 
Step and ramp inputs. 

The results of simulation on an electronic differential analyzer show a 
large reduction in both response time and overshoot for a servo with prediction 
as compared to the same servo with proportional plus error rate control. 

The introduction of an approximation to the true switching curve and 
moderate amounts of friction had little influence on the operation of the pre- 
diction control. The shape of the curve could be altered to optimize for any 
particular combination of Coulomb or viscous friction coefficients; however, the 
improvement was small. 

The application of the nonlinear control to an instrument servo consisting 


iii 





of a 400 cycle per second, two-phase motor coupled directly to a one-turn film 
potentiometer substantiated the results obtained by simulation. The character- 
istics of the signal from the film potentiometer permitted use of an approximate 
differentiator in the control circuit. With this system a resolution of one part 
in 40,000 and a response time of 0.06 seconds to a half-turn step input were 


obtained. 
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SYMBOLS 


Coulomb friction torque, lb. -ft. 

error rate coefficient, 2 //o, 

control voltage 

viscous friction coefficient, lb. -ft. -sec. 
subscript meaning input 

moment of inertia, slug- ft? 


a constant, y 72350 or z /2(Tra + C) 


1 S Ir subscript meaning maximum 


subscript meaning output 
resistance, megohms unless otherwise specified 
time, seconds 
torque, Ib. -ft. 
voltage 
WITZ € 
error, 9; DOS radians 
error to torque-saturate the servo, Tm 
de /dt 
maximum error rate, T,,/f 
de/d? 
damping ratio 
angular position, radians 
gain constant 


3.14159 

















OPTIMUM PREDICTION FOR A 


TORQUE SATURATING SERVO 


INTRODUCTION 


The design of the controller for a servomechanism is based to a large 
extent on the particular type of input function and the requirements on the servo 
such as response time, accuracy, and amount of overshoot. Linear methods of 
control are extensively used as they are easy to implement and analyze. Since 
linear control does not necessarily represent the optimum design in most cases, 
considerable work has been done recently to improve performance by nonlinear 
techniques. 

One such technique involves use of prediction or switching circuits intended 
to give optimum time response and limitation of overshoot for step or ramp in- 
puts or combinations of the two. Previous investigations oí methods of control, 


* 
using some measure of prediction include nonlinear damping functions (50 ; 


prediction with relay servos(C: D) and multiple mode techniques (Ë ) 
Of particular interest is the work of Hopkin(F ) ‚ whose approach to prediction 
is very similar to that in this study although mechanized in a different manner. 


It is the purpose of this investigation to extend previous work by the 


* Parenthetical superscripts refer to references. 








2 
examination of use of a particular type of prediction in the control of a torque 
saturating, second order servomechanism. The basic theory was originally 


developed by L. L. Rauch andR. M. Howe) 


and is summarized in Appendix I. 
This study develops a method of implementation of the theory for simulation in 
several phases and application to a specific instrument servo. 

This investigation was undertaken by William W. Gay and Wayne S. McCord 
as a joint thesis project to satisfy, in part, the requirements for a Master of 


Science Degree in Aeronautical Engineering at the University of Michigan. 


SIMULATION 


Simulation was undertaken to obtain computer solutions to theoretical 
equations of the servo problem under study. These solutions were obtained for 


several variations of the servo setup in order to make a series of comparisons. 


PROCEDURE 


The servo setups considered were those of: 
1. A servo using proportional and error rate control without prediction. 
2. A servo using proportional and error rate control with prediction. 
3. The same servo as in #2 but with an approximate prediction curve. 
a. No friction. 
b. With Coulomb and/or viscous friction 
(but no compensation for friction in the prediction curve). 
c. With Coulomb and/or viscous friction 
(with compensation for Coulomb or viscous friction). 
4. The same servo as in #3, a., but with approximate differentiation of 


error to obtain error rate. 








3 
The first point considered in the simulation was to show the improvement 
in response time and reduction of overshoot by adding prediction to a servo with 
proportional plus error rate control. 
The equation for a torque-saturating, frictionless servo employing proportional 


and error rate control is 
i e 2 に ucl dm P 


To adapt the equation to the computer for solution, it was time scaled by the 


undamped natural frequency, W p? and magnitude scaled by the error required 


for torque saturation, €, . These quantities are defined as 


ZU ge 
ONE 


= 2J 
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Also, C 
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where f is the damping ratio. The equation can then be written 


slo: 


- [e +298] 


sat +6, 


The arbitrary values of parameters chosen were Del f= 0.5, and 

€, = 0.025, where unity was represented by 100 volts in the computer. The 
solutions were then presented in dimensionless form as E/wn versus €/é Í in 
the phase plane and E/E, versuswnt in the time plots. Figure 1 shows the 
computer circuit for solution of this equation. 


The effect of torque saturation was obtained by limiting with operational 


amplifiers. The output stage of each amplifier was a cathode follower. When 
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the outputs of two amplifiers were connected in parallel, the more positive signal 
determined the output of the pair. Thus by introducing a constant positive voltage 
into the input of one of the paralleled amplifiers, the output of the pair was limited 
in the negative direction. By using two such pairs in cascade, both plus and minus 
limiting was obtained. 

The addition of prediction modifies the equation of motion of the above servo to 
the form developed in Appendix I. 

16,= ple + (C8) see2 Ip + な lele| ES 

» 

The method for adapting this equation to a suitable form for computer solution is 


identical to the procedure used previously with the result 


a 
cj, 
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set t€, ee, s satte, 
This equation was actually solved by two different computer circuits. Figure 


2 shows the computer circuit for exact solution. A servo multiplier was used to 


form the product Es € 
So, 





The computer circuit in Figure 3 was also used to solve this equation. In 
this case, a diode limiter network at the input to an operational amplifier generated 
an approximation to the prediction curve. The limiting network and approximate 
prediction curve for this case with no friction are shown in Figure 5. The diode 
networks used in this and subsequent computer circuits were developed by the 
methods presented in Appendix II. 

Friction was introduced into the servo problem with no change in the pre- 
diction curve. 


With Coulomb or viscous friction, the equations, when adapted for computer 








solution, became respectively 
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€), 


Po 
satte, á 
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where T - l/fis the servo time constant, and C is the Coulomb friction torque. 
The effect of Coulomb friction was introduced in the computer as a change 
in torque saturation levels depending on the sign of or Viscous friction was 
introduced in the computer as a feedback in the 6 cQ integrator. The circuit 
additions are shown in Figure 3 in dashed lines. 
Next, compensation for Coulomb and viscous friction effects was introduced 
by changing the approximate prediction function. 
When the prediction is optimized for a specific Coulomb friction, the equation 


for the predictor servo, when adapted for computer solution becomes 


L- -É $i] + 
ae = É +r Bact: Ce EX 


The diode limiting network and approximate prediction curve designed to include 

compensation for Coulomb friction of C/ E = 0.3 were unchanged from Figure 5. 

The only requirement was to multiply the coefficient of € /w_ by the factor : 
y requi ply /o , by Y ITC/TL ラー 


When the prediction curve is optimized for a specific viscous friction, the 


equation, when adapted for computer solution becomes 
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The diode limiting networks and approximate prediction curve compensating for 

viscous friction are shown in Figure 6 for T Wr = 10 and T0) n = 8; 95. 
Approximate differentiation was introduced only for the servo with no friction. 

The computer equation developed previously still applies. The computer setup 


and approximate differentiation circuit are shown in Figure 4. 
RESULTS AND DISCUSSION 


The results of the simulation study were derived entirely from the phase plane 
and time plots for the various conditions outlined previously. Since some effects 
are more apparent in one representation than the other, certain solutions were 
recorded in the phase plane and also as a time function. Complete duplication of 
phase plane and time plots was not attempted. The solutions were grouped so a 
comparison of certain methods of control could be made to determine the answers 
to five questions. 

1. Does the use of this type of prediction materially improve the performance 
of a servo system? 

2. Does an approximation to the prediction curve materially degrade servo 
performance? 


3. How much does friction affect servo performance when prediction is used? 








4. Does compensation for friction effects materially improve the servo 
performance ? 

5. Is approximate differentiation of error to obtain error rate feasible for 
this system? 

A measure of the improvement in performance due to prediction is available 
from the solutions in Figures 7 through 13. For small step inputs the servo with 
prediction shows no reduction in the size of the overshoot or response time. 
However, as the size of the input was increased, the improvement became very 
pronounced. For a step of € / €, s 32, the overshoot was reduced by about 
95% and response time by 65%. 

Of interest in this comparison are the responses, in Figures 14 through 16, 
of the servo with and without prediction to sine wave inputs. Although there is no 
particular significance to this type of response and the prediction was not optimized 
for it, the improvement was still evident. For the highest frequency input, the 
servo with prediction could follow the signal, while the servo without prediction 
developed a divergent error. This condition was shown in Figure 16 where the 
input frequency was 1/3 of the servo natural frequency and the ratio of the maximum 
acceleration of the input to the maximum acceleration of the servo was 8/9. 

Typical solutions of the servo response using an approximation to the pre- 
diction curve are presented in Figures 17 (a) and (b) and 18. A comparison of these 
solutions to those in Figures 9 (a), 8 (a) and 12 respectively shows no preceptible 
difference. Since the operation of the system with this five- segment approximation 
was satisfactory, no solutions were made using a more exact approximation to the 
true switching curve. Allthe remaining solutions were obtained with this method 


of generating the prediction curve. 








The selection of friction coefficients to demonstrate the effects of friction 
was somewhat arbitrary. The Coulomb friction to torque ratio corresponds to 
that specified for the servo later tested. Since the viscous friction of this 
system was small, corresponding to T E 100, additional values of viscous 
friction were selected to emphasize the effect. The change in the response with 
either type of friction is most evident in the phase plane plots of Figures 19 and 
20. This change appears as a considerable reduction in error rate, but it may 
be observed in the error versus time plots in Figures 21 and 22 that the actual 
increase in response time was small even for these large amounts of friction. 

It may be noted that when either type of friction was present the solution path 
lay somewhat within the linear region during deceleration, indicating that less 
than maximum torque was being used. 

In view of the fact that friction has the effect of slowing servo response to 
some extent, it was necessary to determine if response could be improved by 
changing the prediction curve to permit utilization of maximum torque during 
deceleration. 

Figures 23 and 24 show phase plane and time plots where the prediction 
curve was altered to take into account the presence of Coulomb friction. These 
plots may be compared to those in Figures 20 and 22 for the same conditions 
but with prediction calculated for no friction. A very minor improvement in 
response time was obtained by altering the prediction curve to compensate for 
Coulomb friction. More improvement was noted as size of the step was increased. 

In considering the various runs in Figures 25 through 29 it immediately be- 
comes obvious that little improvement in servo time response was obtained when 


the prediction curve was altered to take into account viscous friction. It may be 








noted that Figures 25 and 26 show a case where no velocity saturation due to 
viscous friction occured. Figures 27 through 29 show a case where velocity 
saturation was approached. 

It may be said that prediction to include either type of friction effects is 
feasible but offers little advantage. 

By examining Figure 30 it is apparent that approximate differentiation was 
satisfactory in that the phase plane plot appears quite similar to phase plane 
plots where the true value for error rate was used in computation. 

Stability was a problem in the simulation, however, and would presumably 
be a serious problem in application. The time constants of the differentiating 


circuit for the simulation were selected in order to avoid oscillation. 
APPLICATION TO AN INSTRUMENT SERVOMECHANISM 


A physical system was constructed to demonstrate in practice the use of 
approximate differentiation and an approximate prediction curve similar to the 


one examined in the simulation. 
EQUIPMENT 


Motor. The servo motor was a two phase, 400 cycle per second, 115 volt 
Transicoil motor, type 18M. Rated stall torque - 3. 0 oz. -in., moment of inertia - 
0. 0197 oz. ^in? , stall watts per phase - 18 watts, Coulomb friction - 0.10 oz. -in., 
coefficient of viscous friction - 0. 00285 oz. -in. -sec. 

Potentiometer. Feedback voltage was obtained from a three section, con- 
tinuous film potentiometer, model number 205, made by Computer Instruments 


Corporation, Resistance - 50,000 ohms per section. Total moment of inertia - 
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0. 164 oz. -in. ^ 


, total Coulomb friction - 0.895 oz.-in., total viscous friction - 
negligible. 

Diode limiting circuits. This is the circuit described in Appendix II, 
Figure II-3. 

Compensation and control circuits. The compensation and control circuits 
utilized the operational amplifiers of a standard analog computer. 

Power Amplifier. Two different power amplifiers were used to drive the 
control phase of the motor. One was a115 volt, 400 cycle per second, push- 
pull, magnetic amplifier, Magnetic Amplifiers Inc. Model MA 41501-CY. 

The other amplifier was an A.C. electronic amplifier used in conjunction 
with a 400 cycle per second chopper to convert the D.C. control signal to a 
400 cycle per second A.C. signal. This amplifier was originally designed as a 
60 cycle per second servo amplifier and was adapted for use with this system. 

The choice of equipment was governed primarily by the requirements of the 
control circuits, specifically the differentiator. The film potentiometer was 
selected in preference to a wire wound potentiometer because the signal obtained 
had less noise and better resolution. The wire-to-wire discontinuities of a wire 
wound potentiometer would have introduced large errors in the differentiated 
signal and contributed to the instability of the differentiating amplifier. 

The motor was selected from those available on the basis of largest torque 
to inertia ratio. It was coupled directly to the potentiometer to eliminate the 
backlash and binding problems of gear trains. Although the motor-to-load 
inertia balance was not good, the acceleration loss due to direct coupling was 
small. The best gear ratio, if used, would have been about three to one. 


Initially a magnetic amplifier was chosen for power amplification because 








11 
of simplicity and reliability. The severe bandwidth limitations were alleviated 
somewhat by compensation as described in Appendix H so that the response was 
flat out to 200 cycles per second. However, the response time of the servo was 
short enough that the 200 cycles per second bandwidth still represented a limi- 
tation which appeared as a time delay. When the system gain was increased to 
obtain reasonable resolution, the time lag in the magnetic amplifier sustained 
oscillations of the servo at a frequency of 80 to 120 cycles per second, depending 
upon the particular gain setting. The oscillations were similar to those that 
would be experienced by an identical bang-bang servo with a fixed switching 
delay of about 5 milliseconds. This time delay was verified by observation of 
the control and output signals of the magnetic amplifier on a dual beam oscillo- 
scope. The electronic amplifier and chopper were substituted when operation 
with the magnetic amplifier did not produce the potential resolution or tracking 


performance of the servo. 


PROCEDURE 


The servo was operated to obtain a comparison of the system with pro- 
portional control, proportional plus error rate control and proportional plus 
error rate control with prediction. The comparison was made by observing 
only the step response, and no attempt was made to analyze the servo completely. 
AS an item of interest, the velocity lag and constant input errors were measured. 

The basic servo circuit diagram is shown in Figure 31. Appropriate modi- 
fications to this basic circuit were made to obtain the particular types of control. 
The servo loop gain was adjusted to the highest possible value consistent with 


stable operation. If the gain was further increased, the servo went into a very 
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low amplitude oscillation at about 160 cycles per second. 

The time constants of the approximate differentiating circuit were determined 
empirically with a view toward limiting noise and high frequency oscillations yet 
retaining sufficient bandwidth for good operation. The time constant of 0.0005 
seconds corresponds to 319 cycles or roughly twice the maximum response 
frequency of the servo. The basic differentiating circuit, with constant gain at 
high frequencies (within the capabilities of the operational amplifier), passed 
too much noise, so the additional capacitor was added in the feedback to give a 
-6 db per octave roll off of gain for frequencies above 319 cycles per second. 

The approximations in differentiation modified the error rate signal such 
that the amplitude was not increased in proportion to the frequency as prescribed 
for exact differentiation. Some compensation for the resulting loss of gain was 
made in the next amplification stage where the product Valor E was 

(Ti + C) 
formed. In place of using a calculated gain of 3. 16 in this amplifier, the gain 
was increased to 4. The value of 4 was determined empirically for satisfactory 
servo response. The factor T in ars is the scale factor between the 
output measured in radians and computer units, where 100 volts represents unity 
in the computer. 

Photographs of the cathode ray oscilloscope presentation were made for 
phase plane and error versus time response of the servo to a 100 volt peak-to- 
peak 3 cycle per second square wave input. Similar photographs were made of 
the error versus time response for a 10 volt peak-to-peak 5 cycle per second 
square wave input. 

The constant input error was measured by disturbing the servo slightly 


and noting the maximum error when the output had come to rest. The error 








13 
was actually measured on a voltmeter having a full scale deflection of 0.1 volt. 
The velocity lag was measured in a manner similar to the constant input 
error. The effect of a ramp input was produced by a 0.5 cycle per second 
triangular wave. The frequency of the triangular wave was low enough that the 
transient effects damped sufficiently to permit accurate measurement of the 


velocity lag. 
RESULTS AND DISCUSSION 


Qualitatively the performance of the servo corresponded closely to results 
of the simulation. Although specific performance of the servo was not an 
objective, the response time, ability to track a slowly varying signal and the 
resolution approached closely the potential capability of the system. 

The measured value of maximum steady state error was + 0.005 volts. 

The terminal voltages on the potentiometer were + 100 volts with the center- 
tap grounded. Based on this reference voltage range of 200 volts, the maximum 
steady state error corresponds to resolution of at least one part in 40, 000. 

From the maximum steady state error, an approximation to the width of 
the linear region (2 €,) may be calculated. The Coulomb friction to rated 
torque ratio was 0.33. However, actual measurement of the stall torque of the 
motor as driven by the electronic amplifier was about 2 0z. -in., giving an actual 
friction to torque ratio of about 0.5. Using this ratio, the linear region was 
calculated to be in the order of 0. 02 volts. 

Among the factors contributing to the low value of torque were poor output 
waveform from the electronic amplifier and somewhat less than 90° phase shift 


between the two motor phase voltages. 
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The photographs of the oscilloscope presentations of the phase plane and 
error versus time responses are shown in Figure 32. The system using propor- 
tional control exhibited a series of fairly large overshoots and long response 
time. Addition of error rate control ( Y = 0.5) reduced the size and number of 
overshoots and response time. Furthermore, addition of prediction essentially 
eliminated overshoot and very substantially reduced response time. As expected, 
the improvement due to prediction was much more pronounced for the larger 
Steps. As an example, for the 100 volt step input the response times were 0.17, 
0.13 and 0.06 seconds for the three types of control. 


Measurements were made of velocity lag and are recorded in Table I. 


TABLE I 


VELOCITY LAG (All Types of Control) 


01 | 

014 
.02 
03 
.04 
.06 
.076 
. 08 
.10 
13 












どど どど と ど と どど どら どど の どの 


Since the rotation of the motor was limited to at most one revolution, the 
velocity never exceeded 20% of the maximum motor velocity. Within this range, 
the viscous friction effects of the servo were negligible. For practical con- 
siderations, this fact plus the very narrow linear region makes the response of 


the servo to a step input of several volts very similar to that of a bang-bang servo 








> 
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with only Coulomb damping. 

The response of the servo with proportional control to the 100 volt step 
input exhibits a ratio of successive overshoots of about 0.4. Calculations 
using this ratio for a theoretical bang-bang servo with Coulomb damping 
indicate a Coulomb friction to torque ratio of 0.43, which checks closely with 
the ratio of the measured values, 0.5. 

The response of the servo with prediction can also be compared to the 
theoretical response of a corresponding bang-bang servo. Again for the 100 
volt step input, the observed response time was 0.06 seconds which compares 


to a computed time of 0. 063 seconds. 


CONCLUSIONS 


Based on the results obtained in the simulation and in application to the 

instrument servo it is concluded that: 

1. Use of prediction reduces response time and overshoot, particularly 
for large step and ramp inputs. 

2. Use of prediction improves transient response of the servo for arbitrary 
input signals, such as a sine wave input, even though the prediction was 
optimized only for step and ramp inputs. 

3. Use of simple approximations to the exact prediction curve does not 
degrade performance. 

4. Modification of the prediction curve to optimize for friction effects 
was possible but little improvement was achieved over the prediction 
curve optimized for no friction. 


5. Noise level in the output signal from the film potentiometer was low, 








16 
thus permitting differentiation with an operational amplifier to obtain 
approximate error rate signal. 

Resolution oí the film potentiometer and direct coupling permitted a 
degree of control sufficient to give servo resolution of one part in 
40,000. 

For high loop gains the region of linear operation was very narrow, 
hence, the system acted much like a bang-bang servo with Coulomb 


damping. 
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FIGURE 5 — (a) Approximate Prediction Curve for a Frictionless 
Servo, K = M/2T m» and a Servo with Coulomb 
Friction, K = P/UT mn + C) 
(b) Corresponding Diode Function Generator Circuit 
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FIGURE 6 — (a) Approximate Prediction Curve for a Servo with 
Viscous Friction, 79 — 10, 8.95 
(b) Corresponding Diode Function Generator Circuit 
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IGURE 7 — Step Response of a Frictionless Servo with 
Proportional and Error Rate Control, C, = 1/w, 
(a) With Prediction (b) Without Prediction 




















FIGURE 8 — Step Response of a Frictionless Servo with Proportion: 
and Error Rate Control, C, = 1/o) 
(a) With Prediction (b) Without Predictior 
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FIGURE 9 — Step Response of a Frictionless Servo with Proportional 
and Error Rate Control, COS 1/0, 


(a) With Prediction (b) Without Prediction 
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FIGURE 10 — Step Response of a Frictionless Servo with 
Proportional and Error Rate Control, C, = lon 
(a) With Prediction (b) Without Prediction 
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FIGURE 11 — Ramp Response of a Frictionless Servo with 
Proportional and Error Rate Control, C, - 1/w n, 
9 i/ 9g = 0.1 
(a) With Prediction (b) Without Prediction 
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FIGURE 14 — Sine Wave Response of a Frictionless Servo with 
Proportional and Error Rate Control, Ce = 1⁄2n, 
" - 1/5, €, (max)/ Ó, (max) - 0.32 
(a) With Prediction (b) Without Prediction 
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FIGURE 15 — Sine Wave Response of a Frictionless Servo with 
Proportional and Error Rate Control, C¿ = 1/0 p, 
42/00 =1/4, 0, (max)/ O (max) = 1/2 
(a) With Prediction (b) Without Prediction 











FIGURE 16 — Sine Wave Response of a Frictionless Servo with 
Proportional and Error Rate Control, C¿ = 1/0 y 
a = 173, Ó, (max) / 0, (max) - 8/9 
(a) With Prediction (b) Without Prediction 
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FIGURE 17 — Step Responses oí a Frictionless Servo with 
Proportional and Error Rate Control, C¿ = 1/0 p, 
Plus Approximate Prediction. 
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FIGURE 21 — Step Response of a Servo with Prediction Optimized 
for No Friction 
(a) No Friction (b) f»4 2100 (c) TOR = 10 
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FIGURE 22 — Step Response of a Servo with Prediction Optimized 
for No Friction (a) No Friction (b) No Viscous 
Friction, C/T,, = 0.3 (c)7«w p = 10, CARS 
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FIGURE 24 — Step Response of a Servo with Prediction 
Optimized for Coulomb Friction, C/T n = 0.3 
(a) No Viscous Friction (b) To p = 10 
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FIGURE 26 — Step Responses of a Servo, Tw  - 10, with Prediction 
: n 
Optimized for 


(a) No Friction (b) Viscous Friction, UM = 10 
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FIGURE 29 — Step Responses of a Servo Approaching Velocity 
Saturation, Ta à = 8.95, Prediction Optimized for 


(a) No Friction (b) Viscous Friction, T% n = 8-95 
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APPENDIX I 


OPTIMUM SWITCHING CRITERIA 


The minimum response time for a servo system with limited torque is 
obtained if the maximum available torque is applied to accelerate the system, 
and then the torque reversed for maximum deceleration to arrive simultaneously 
at zero error and zero error rate. If the torque is then removed, the error will 
remain zero until another command signal is given to the system. This type of 
response is discussed for three categories; a system with no friction, a system 
with only Coulomb friction, and a system with only viscous friction. 

The problem of determining the exact time to reverse the torque depends 
in part upon the characteristics of the command signal. For the first two cate- 
gories the acceleration of the command signal is assumed to be zero. That is, 
the switching time is optimized for constant or uniformly varying command 
signals. For the case with viscous friction, the switching time is optimized 
for a constant command signal. 


First consider the system with no friction. The equation of motion is 
I o5 z+ Tia = torque 


or, in terms oí error, since 6, = 0, 


Tm 
I 


QN 





Integrating two times and assuming initial conditions of E(0) = 0 and 
€ (0) = 0, the solution of the equation is 
2 


T t 
I I 





€ = + 








1-2 


or 





since t - 4 1 € from the first integration. 


m 
In the phase plane of error and error rate, the solution paths are two 


parabolas passing through the origin as shown in Figure I-1. 


Tm < O 





FIGURE I-1 


Phase Plane Plot of Solutions Through the Origin 


The parts of interest lie in the second and fourth quadrants, the other two portions 
represent divergent solutions. The phase plane is thus divided into two regions 
by the two segments of parabolas comprising the switching curve. The equation 


for this curve is 


I TEE 
aa = lele 
2T m 











I-3 
For a solution starting above and to the right of the switching curve, full positive 
torque applied will cause the solution to move along a parabolic path until it 
reaches the switching curve in the fourth quadrant. If the torque is then reversed, 
the response will follow the switching curve to the origin, and the torque can be 
reduced to zero. The same procedure applies to the other half of the phase plane 
with the opposite sense of torque. Sample phase plane trajectories are shown in 


Figure 1-2. 





FIGURE 1-2 


Example Trajectories 


Although this switching criterion represents the fastest response for large 
errors or error rates, a practical system mechanized in this manner would not 
follow slowly varying inputs smoothly due to switching delays and other non-ideal 
characteristics. A linear system with error rate or derivative feedback in 


addition to proportional control is better for this type of operation near the origin 











I-4 
in the phase plane. The torque in a linear system with proportional plus error 
rate control is proportional to € +C, é , and the optimum nonlinear system 
utilizes maximum torque tending to satisfy the equation 

c jo er 
2 Tm 
In combining the two systems, the torque should be made proportional to 
E + C。6 ae — . Near the origin of the phase plane, € is small 
and the term _I lé|é is negligible compared to the others, hence the servo 


E 
response is linear as desired. Away from the origin the C¿€ term must be 





limited so the response will tend to make 
C + ラキ ー I e z 0 
Tm 


For a large gain constant, LI a small error, € z エル バン , will torque 

saturate the system. If the C。 と term is made to saturate at a level exactly 
equal to the error required to torque saturate the servo, then away from the origin 
the response will lie exactly on the optimum switching curve. The equation of 


motion of the servo is then 


I6,= u|e + (Ca Ekar: tar E 
er 2 sat: 2 


where sat + Tan A represents the limits on the particular quantity in brackets. 

In the phase plane, the torque equation represents a switching curve with a 
linear region adjacent to the switching curve. The limits of this linear region 
are shown in Figure 1-3. The zero torque line lies midway between the boundaries 


of the linear region in the no friction case. 











FIGURE I-3 


Phase Plane Trajectory For the Optimum 


Switching Curve with Linear Region 


The extension of the above development to the case of a system with Coulomb 
friction is quite straightforward. The Coulomb friction represents a constant 
torque, C, which adds to servo torque during deceleration and subtracts from 


servo torque during acceleration. The switching curve equation then becomes 


I 


E + anr O 


lé | € =. 0 
and the corresponding equation of motion of the system is 


.. ñ E I “le 
E Dx 


X. 
sat t =m 
S 


The effects of Coulomb friction on the switching curve are to increase the 


magnitude of the slope in a ratio (1+ C/T,,), and to shift the zero torque line 








from the center of the linear region. 


In developing the equation of the switching curve for a system with viscous 
friction, the additional condition of 9 i= 0 was applied. The equation of motion 


of the system is now 
Ic f 66 = torque = = Ta 
or in terms of error 
Te + と = 37€ 
m 
where T - l/f- time constant and c T, /f = maximum error rate 


ŻE re | , (Q) = 0 


or t= -Tln(1£ é / É m) 


im oe 


Ez :zég|t-Tü-e €(0) - 0 


E= t ém T a EE EN 
€rré =-T é, mü + | é /¿ ime 656 


€+Té= -Té 1 + |é /é |), €<o 


The equation of motion of the servo system is then 


Iö,+r6, = “le *(C, €) 


sata Je ji cis on EL m 





1-7 
To compare the switching curve for viscous friction to the one with no 


friction, consider for the moment, positive é. 
Letx = 1 T/ Z Ta E 


The switching curve for no friction is 


The switching curve for viscous friction is 


E -Té,| ina« €/€..) - eye] 


と = 2/m* | ma + mx) - mx | 
Since mx <1, expand In(1 + mx) in a series 


と = 27[ [mx- al, tmx”. mar, ih m te .] - mx) 


= 2 3 4 
2 (mx) (mx) (mx) 
ezda | 52 + S - — + | 
2 3 
ao +e | 


It is apparent that for no friction, m = 0, the two curves are identical. 
The difference is important only when mx is significant compared to 1. The same 


result can be obtained by comparing the slopes of the two curves. 


SE- -2x forf - 0 


and 
de_ -2x 
& = T+ mx forigo 





Comparative curves are shown in Figure I-4 for several values of m. 


The curves are terminated at mx < 1. 
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APPENDIX H 


THE FUNCTION GENERATOR 


Any single valued function can be approximated by straight line segments 
to a reasonable degree of accuracy with an operational amplifier and biased 
diodes. The problem is only one of calculating the resistance values and bias 
voltages required to give the proper slope and location of the straight line 
segments. 

Consider the circuit and response curve in Figure II-1. Make the assumption 
that the resistances are large compared to the grid to plate resistance of the 
diode and the power supply impedance so that the voltage drop across the con- 


ducting diode and in the power supply are negligible. 





TEST Rit Bs 
O - Ry v) ej = v 
R, + R 
es h] ei ei ご - + v 
o = 1 
R; + Ro 
FIGURE II-1 
(a) Diode Limiting Circuit (b) Response Curve 


If another input is now added with the diode reversed and biased positively, 


the output can be made symmetrical with respect to the origin. Symmetry is 








H-2 
not necessary, but is demonstrated below in Figure II-2 because of the par- 


ticular function of interest in this report. 


+ 


i A 





1 2 
es = 0 leil < = v 
R1 t Ro 
nu e Jeil > R. v 
1 
FIGURE H-2 
(a) Symmetrical Diode (b) Response Curve 


Limiting Circuit 


By adding other inputs with different resistances and bias voltages, more 
straight line segments can be generated and combined with the existing ones. 
Each input contributes a change in slope in the region where the diode is cut 
off, and a constant output elsewhere. Thus, by choosing the proper cut off 
voltage and resistance, the point of slope change and the amount of change 
are determined. 

For most problems it is desirable to have a well defined, steady zero 
reference. If the diodes are conducting for zero input, drifting of the bias 


voltage, changes in the diode resistance and inequalities of resistances tend 





H-3 
to give output drift and balancing problems. For these reasons it is preferable 
to operate with the diodes cut off at zero input. This condition in general re- 
quires an additional operational amplifier if the second derivative of the function 
is positive, as it is in this servo problem. The procedure will be demonstrated 
for generation of the torque switching curve used for the frictionless servo. 

It should be noted at this point that biased diodes could be used in the 
amplifier feedback circuit with an effect opposite to those in the input. This 
procedure was not used because it is difficult to generate a zero slope, and also 
a change of slope of any segment involves solution of a set of "n" simultaneous 
equations to determine the new resistance values where (2n - 1) is the number 
of segments in the curve. 

The mechanization of the switching curve involved generation of the functions, 
ーー | ë | com for the application to the servo, and the equivalent, 


E | <_| € _, for the simulation. In the first case € was multiplied by 
2Tm Wy Da 


1 V27Tw  andinthe second case € / wp Was multiplied by 4 AT : 
The multiplication conveniently scaled the input, e, and also reduced the problem 


to one of generating |e,| e; for either case. The curve was approximated by five 
segments from -100 Kei < 100 volts. This range corresponds to that from minus 
one to plus one in computer units. More segments were tried, but the change in 
performance of the servo was not noticeable. 

The steepest slope was generated by one amplifier and the output of the other 
amplifier subtracted from it to generate the desired curve. The slopes used, 
range of input voltages, resistance, and bias voltages are tabulated below for 
positive e; . The negative half of the curve is symmetric with respect to the 


origin. A one megohm feedback resistor is used. 





H-4 
TABLE H-1 


Diode Function Parameters (No Friction) 





ei(volts) Slope | Change of Slope Ry, = Rg Bias(Volts) . 


The resulting circuit appears in Figure II-3. 


tov. -5u% 





il 


+ POW -pOu 


FIGURE 1-3 


Diode Function Generator (No Friction) 


The outputs of the two amplifiers are shown below in Figure II-4 along 


with the contribution of the -ei input to amplifier B. 








FIGURE II-4 


Diode Function (No Friction) 


Altbough this particular approximation appears very crude, a comparison 
of it to the true function shows there is no more than one volt error over half 
the range and four volts over 90% of the range. 

By adding one more segment on each end, the error can be kept within one 


volt over 97. 5% of full scale. 





APPENDIX III 


MAGNETIC AMPLIFIER CHARACTERISTICS AND COMPENSATION 


To determine the static characteristics of the magnetic amplifier, a 
control current was supplied from an operational amplifier and the output 
terminals were connected to a 1320 ohm resistance load. This load represented 
ten watts at rated voltage. The static response curve is shown in Figure III-2. 

The frequency response of the magnetic amplifier was obtained by again 
driving it from an operational amplifier with current feedback. The operational 
amplifier offered no bandwidth restrictions over the range of frequencies con- 
cerned. The frequency response curve in Figure III-3 shows that the output was 
down 3 db at about 17 cycles per second. By adding a 1.6 uf condenser and 
900 ohm resistor in parallel with the 10k resistor to ground, the output was 
kept quite flat out to about 200 cycles per second after which it falls off again. 


The control circuits are shown in Figure III-1. 


(a) 





FIGURE III-1 


Magnetic Amplifier Control Circuits 


(a) Uncompensated (b) Compensated 
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